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The phonon dispersion was measured at room temperature along (0,0, L) in the tetragonal phase 
of LaFeAsO using inelastic x-ray scattering. Spin-polarized first-principles calculations imposing 
various types of antiferromagnetic order are in better agreement with the experimental results than 
nonmagnetic calculations, although the measurements were made well above the magnetic ordering 
temperature, Tn- Splitting observed between two A\ g phonon modes at 22 and 26 meV is only 
observed in spin-polarized calculations. Magneto-structural effects similar to those observed in the 
AFe2As2 materials are confirmed present in LaFeAsO. The presence of Fe-spin is necessary to find 
reasonable agreement of the calculations with the measured spectrum well above Tn - On-site Fe 
and As force constants show significant softening compared to nonmagnetic calculations, however an 
investigation of the real-space force constants associates the magnetoelastic coupling with a complex 
renormalization instead of softening of a specific pairwise force. 



PACS numbers: 74.25.Kc, 78.70.Nx, 74. 25. Ha 

Despite rather convincing arguments that supercon- 
ductivity in the AFe 2 As 2 (A=Ca,Sr,Ba,Eu) and RFeAsO 
(R=La,Ce,Pr,Nd,Sm,Gd)-based compounds does not 
originate from conventional electron-phonon coupling, 1 
these systems do display significant sensitivity to the lat- 
tice geometry. For example, the size of the Fe moment 
is sensitive to the lattice parameters and As position, as 
shown by Density functional theory (DFT) calculations. 
One thus expects strong magneto-structural coupling in 
these compounds!^ Also, measurements of the room tem- 
perature phonon density-of-states (DOS) in LaFeAsO in- 
dicated some disagreement with non-spin-polarized DFT 
calculations Distinct features in the phonon DOS, 
likely associated with atomic displacements in the Fe-As 
plane, were observed at significantly lower energies than 
non- magnetic calculations suggest. It was noted (em- 
pirically) that softening of the Fe-As force constants by 
30% brings the calculated phonon DOS into better agree- 
ment with the dataP Theoretical studies have shown that 
strong coupling between Fe magnetism and the As posi- 
tion leads to the softening of the Fe-As force constants, 
thereby explaining the observed phonon spectraP 

While these magnetostructural effects are well docu- 
mented in the AFe2As2-based systems, it is not clear 
if the same effects are present in the RFeAsO system. 
One key example of this coupling in CaFe2As2 comes 
from the observation of a transition from the antiferro- 
magnetic state to a non-magnetic "collapsed tetragonal" 
state under applied pressure. 7 In this case, a reduction of 
the c-axis lattice parameter by 9.5% is associated with 
the complete collapse of the Fe magnetic moment P 

The lattice vibrational frequencies associated with c- 



axis vibrations of Ca and As atoms in CaFe2As2 and 
BaFe2As2^^ have been shown by inelastic neutron and 
x-ray scattering to disagree with predictions of non-spin 
polarized DFT calculations. In particular, the energy 
splitting between c-axis phonon branches containing As 
displacements was found to be in strong disagreement 
with non-spin polarized calculations. Ultimately, spin- 
polarized calculations in the local spin density approxi- 
mation that include the AFM order present at lower tem- 
peratures were required to bring the calculated phonon 
dispersion into better agreement with room tempera- 
ture measurementsP Our group was able to confirm 
the role of magnetism in c-axis polarized Ca and As 
modes in CaFe2As2 using single-crystal inelastic x-ray 
scattering (IXS) measurements at the Advanced Pho- 
ton Source (APS) in combination with spin-polarized 
calculations using the Perdew-Burke-Ernzerhof (PBE) 
exchange-correlation functional.^ 

It might be expected that the presence of RO lay- 
ers, which results in a larger spacing of the FeAs layers 
along the c-axis, might mitigate these effects to some de- 
gree. However, the difficulty of synthesizing RFeAsO in 
single-crystalline form has prevented a quantitative con- 
firmation of similar magneto-structural coupling across 
the AFe2As2 and RFeAsO systems. Recently single crys- 
tal samples of LaFeAsO have become available.^ IXS 
phonon data was recently reported on PrFeAsOo.sP and 
SmFeAsO 13 single-crystals, however the dispersions were 
only measured along the (100) direction. Also, the role 
of spin-phonon coupling could not be ascertained since 
results were compared only to non-spin-polarized LDA 
calculations. 
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TABLE I. Theoretically relaxed and experimentally observed 
z-position for La and As atoms, and the associated magnetic 
moment per Fe atom and total energy. In each case the exper- 
imental l attice parameters of (a = 4.03533A, c = 8.74090A) 
were used™ 



ZLa 
ZAs 

E(Ry) 



NM SDW Striped Checkerboard ExpP^ 

0.13993 0.13875 0.13883 0.13887 0.14154 

0.63829 0.64820 0.64770 0.64401 0.6512 

0.0 2.32 2.30 1.91 0.36-0.78 

0.0 -0.032 -0.033 -0.009 



LaFeAsO single crystals were synthesized in an NaAs 
flux at ambient pressure as described elsewhere.^ Inelas- 
tic x-ray scattering measurements were performed on the 
HERIX instrument at sector 30-ID-C of the Advanced 
Photon Source at Argonne National Laboratory with in- 
cident beam energy of 2 3.72 4 keV and with an energy 
resolution of 1.44 meVP^^ Scattering is described in 
terms of the tetragonal P4/nmm unit cell where Q = 
^ {hi + fcj) + ^f'k. The vectors i, j, and k are the 
fundamental translation unit vectors in real space. Be- 
low Ts=156K, the sample transforms to a n orthorhom- 
bic structure with space group Cmma ! 16 * 17 ^ The relation- 
ship between the Miller indices in the tetragonal P4/nmm 
and orthorhombic Cmma phase are, h = (H Q + K Q ) , k — 
(H — K ), and I = L Q . Below the magnetic ordering 
temperature T/v=138K, the sample develops long-range 
spin-density wave (SDW) AFM order. The sample was 
mounted in the (hhl) plane in a displex for low temper- 
ature studies, and the displex was attached to a 4-circle 
diffractometer. 

Based on previous studies of c-axis polarized phonons 
in CaFe2As2, we focused our study on phonon branches 
along the (0, 0, 8 + £) direction in the Brillouin zone. In 
order to study the dispersion and potential line broaden- 
ing of the phonon modes, the scans were fit to several 
peaks using a pseudo-Voigt line profile. The normal- 
ized pseudo-Voigt function is given in Eqn. [I] where 
fa (x\ T) and /l (x\ V) are normalized Gaussian and 
Lorentzian functions respectively. The mixing parame- 
ter r\ = 0.74, and resolution full- width- at-half-maximum 
(FWHM) T = 1.44 meV was determined from fits to the 
elastic scattering width of Plexiglas. 



f pV = (l- V )f G (x;T)+ V f L (x;T) 



(1) 



Figure [T] shows a line scan consisting of several phonon 
excitations at Q = (0,0,8.3) and Q = (0,0,8.5) at room 
temperature. The peak positions for these and other 
scans were obtained from fits and used to construct the 
dispersion of phonon branches along the different scan di- 
rections, as shown in Fig. [2] The intensity of the phonon 
modes is also represented in Fig. [2] by the diameter of 
the circles. 

In order to understand the features of the phonon dis- 
persion, the experimental measurements were compared 
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FIG. 1. (color online) Energy scan at constant-Q at a) 
Q = (0.0,0,8.3) and b) Q = (0.0,0,8.5) measured at room 
temperature on LaFeAsO. Experimental data are given by 
solid green points. The black line is fit using a pseudo-Voigt 
function. 



to ab initio calculations of the phonons. The phonon 
dispersion was calculated using DFT and Density Func- 
tional Perturbation Theory (DFPT)P There are sig- 
nificant differences in the experimental lattice parame- 
ters and parameters from the "relaxed" structure with 
the lowest calculated energy. Also, in spin-polarized 
calculations with the experimentally observed AFM or- 
der, the lattice distorts into the orthorhombic Cmma 
structure observed experimentally at lower temperatures. 
With these difficulties in mind, the experimental lat- 
tice parameters at room temperature in the tetragonal 
phase (a = 4.03533A , c = 8.74090A) were used for 
all calculations.^^ In addition, there is debate over the 
appropriate internal ^-parameter to use f or the position 
of lanthanum and arsenic atoms lH £l 21 l 22 l For better ac- 
curacy of the calculated phonons, we chose the calcu- 
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lated relaxed positions where all forces were zero. Struc- 
tural parameters used for the non-magnetic and spin- 
polarized calculations as well as experimental measure- 
ments are given in table [IJ The pseudopotentials cho- 
sen used the Perdew-Burke-Ernzerhof (PBE) exchange 
correlation functional ! 23 * 24 * Settings of an 8x8x4 (non- 
magnetic), 4x4x4 (striped & SDW) and 8x8x2 (checker- 
board) k-mesh and 50 Ry and 660 Ry energy cutoffs for 
the wavefunctions and charge density were chosen to en- 
sure meV precision of the calculated phonon dispersion. 
These parameters are s imilar to other phonon calcula- 
tions of LaFeAsO! 25 * 26 ! Phonon frequencies were calcu- 
lated on either a 4x4x2 (nonmagnetic), 2x2x2 (striped & 
SDW) or 4x4x1 (striped) q-mesh and then interpolated 
along several symmetry directions. The resulting phonon 
frequencies and eigenvectors were used to calculate the 
dynamical structure factor along the selected scan direc- 
tions. The dynamical structure factor, which is propor- 
tional to the x-ray scattering intensity, is given in Eqn. 
(2pH^ In these equations, W d (Q) is the Debye- Waller 
factor, rij (q) is the Bose-Einstein distribution, fd (Q) is 
the x-ray form factor, and cr J d (q) is the eigenvector cor- 
responding to the normalized motion of atom d in the 
j th phonon branch. While the DFT calculation does not 
include temperature dependence, the Bose-Einstein dis- 
tribution was set to 300K. 



\H j (Q)| 



H k (Q) = E ex P (~ Wd (Q) + • d) {q • oJ (q) 



M d 



W d (Q) 



h 



4M d Q BZ J a 



(2b) 



Q • cr dj | 



(2 nj (q) + 1) 



(2c) 



Both the x-ray form factor and the Debye- Waller factor 
decrease intensity of the phonon excitations with increas- 
ing Q. The preferred approach for computer applications 
are numerical approximations to the x-ray form factor. 
The x-ray form factor has been parameterized by Waaz- 
maier and Kirfel as the sum of five Gaussians plus a con- 
stant termP^ The Debye- Waller factor, calculated using 
Eqn. |2cJ can be thought of as the mean-squared value of 
the displacement each atom dotted with Q. The vo lume 
integral was calculated using the tetrahedron methocP^^ 
on a 16x16x8 (nonmagnetic), 16x16x4 (checkerboard) 
and 12x12x12 (SDW & striped) Monkhorst-Pack q-point 
gridP^To second order, the integral over the tetrahedron 
is simply the function evaluated at the center point mul- 
tiplied by the volume. To avoid repeating the calculation 
for each value of Q, the nine potential components of the 
phonon eigenvector were stored and the dot product cal- 
culated later. 

The delta function in uo was convoluted with the elas- 
tic scattering width of Plexiglas measured on analyzer 



5. The pseudo-Voigt function fits well to the center of 
each peak, but small discrepancies exist in the tail. To 
minimize this effect a discrete linear convolution between 
the raw experimental data and simulated delta functions 
(single point on a grid) was performed numerically. 

In addition to the energy resolution, the diameter of 
the analyzer leads to a finite resolution in Q. Slightly 
different positions on the analyzer can be described by 
a radial component, determined from the size of the an- 
alyzer (10 cm) and the distance from the sample to the 
analyzer (9 m), and an angular component covering the 
entire circle. Values of Q accepted by the analyzer can 
be written as a function of these two variables. 5000 
samples from a pseudo-random number generator gave 
sufficient precision for convolution of constant-Q scans 
with the Q-space resolution. Due to the lower required 
precision 1000 samples were used for each contour plot. 

Figure [I] shows a constant-Q energy scan at (0,0,8.3) 
and (0, 0, 8.5) at room temperature. Experimental data 
are given by the green points and pseudo-Voigt fits by 
the solid black line. The default values for r] and T only 
account for the energy resolution. At (0,0,8.3) the fit 
on the acoustic mode at 5 meV, however, is adjusted by 
including r] and T as variables. The fitted values of r] and 
T are 0.48 ± 0.07 and 3.96 ± 0.09, respectively. Optical 
modes are present at 22, 27, and 34 meV. At (0,0,8.5) the 
fit on the acoustic mode at 8 meV, however, was adjusted 
by setting r] = 1.0 (Lorentzian function) and including T 
as a variable. The fitted value of T is 3.67 ± 0.12 meV. 
A nearby optical mode is present at 11 meV, along with 
^three other modes at 22, 26 and 34 meV, respectively. 
J Fig. [3] shows several calculations of the dynamical 
structure factor at Q=(0, 0,8.5) which can be directly 
compared to Fig. ^p. The red dotted line is a non- 
magnetic calculation. Frequencies for the acoustic and 
lowest optical modes are reasonable, but the calculated 
intensity of the optical mode is too high. Attempts to 
include the experimental uncertainty in Q could not re- 
produce the observed broadening of the acoustic mode, 
meaning it is not an artifact of Q-space resolution. The 
phonon excitation near 24 meV consists of two modes 
separated by 0.2 meV. At the zone boundary, lower and 
upper Ai g modes consist of As and La motion, respec- 
tively, polarized along the c-axis. At the zone center, 
these modes are mixed, each containing both La and As 
motion and the upper mode contributes 80% of the struc- 
ture factor. This result from the nonmagnetic calculation 
is inconsistent with the measurements, where these two 
modes are clearly split by 4 meV at (0,0,8.5). While the 
calculated frequencies agree with other published phonon 
dispersions j 25 * 26 * they are a few meV lower in energy than 
observed. Small changes in lattice parameters are not re- 
sponsible, as an unphysical 7% reduction in the unit cell 
volume is required to stiffen this phonon mode in the 
nonmagnetic calculation to the observed value. While 
this discrepancy exposes limits on the accuracy of these 
DFT calculations, this should not detract from qualita- 
tive changes between calculations, such as the splitting of 
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FIG. 2. (color online) Contour plots of the calculated dynam- 
ical structure factor along (0,0,L). Values range from blue (no 
intensity) to red (high intensity), and have been multiplied 
by the energy to improve visibility of the optical modes. The 
white dots show the experimentally determined frequencies, 
as described in the text, with the intensity shown by the size 
of the dot. a) Nonmagnetic calculation b) SP calculation with 
SDW ordering, c) SP calculation with checkerboard ordering, 
d) SP calculation with striped ordering. 



the Ai g branches, that are also observed experimentally. 
At both values of Q, the 32 meV feature consists of both 
Fe and As motion, but the intensity is extremely weak. 

In the spin-polarized calculation corresponding to the 
observed stripe AFM structure (Fig. [I^d), the effect of 
the magnetization on Fe is to strongly split these two 
branches at (0,0,8.5) with the 21 meV excitation, con- 
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FIG. 3. (color online) Dynamical structure factor calcula- 
tion of constant-Q line scan at Q = (0,0,8.5). The dotted 
red line corresponds to non-magnetic calculations of the dy- 
namical structure factor. The solid green line corresponds to 
spin-polarized calculations imposing the SDW AFM ordering 
observed at lower temperatures. The black dashed line and 
blue dashed-dotted lines correspond to spin-polarized calcula- 
tions with a striped (ferromagnetic along c) and checkerboard 
ordering, respectively. The experimentally observed frequen- 
cies in Fig. ^3 are shown with vertical grey lines. 



taining As motion, lowering in energy by approximately 
8.6%. The ratio of intensities between the acoustic and 
nearby optical mode moves in the direction of, though 
slightly more than, what is observed experimentally. The 
24 meV peak is primarily La motion. The intensity of the 
32meV feature is 5.2 times stronger, in better agreement 
with experiment. 

In order to better understand the importance of the 
specific magnetic order and size of the Fe moment on the 
lattice dynamics, two additional calculations were per- 
formed in hypothetical magnetic structures. First is the 
"checkerboard" magnetic structure, shown in Fig. |4]3. 
It is a tetragonal space group, where Fe neighbors have 
opposite spins. This calculation converges to a solution 
0.023 Ry higher in energy and energy and with an 18% 
smaller magnetic moment per Fe atom. The acoustic 
mode is slightly softer and has greater intensity. The 21 
meV excitation, containing As motion, is lower in energy 
by approximately 19.9%. The intensity of the 32meV 
feature is 2.9 times stronger than in the nonmagnetic 
calculation. 

Second is the CeFeAsO structure,^ also referred to as 
"striped," and shown in Fig. [4]i. It is an orthorhombic 
space group with ferromagnetic coupling of Fe moments 
along the c-axis. The dynamical structure factor for this 
material is shown with black dashes in Fig. [3j The fre- 
quency and intensity of the acoustic and optical modes 
at 8 and 11 meV are nearly identical. Once again, the 
effect of the magnetization on Fe is to strongly split these 
branches, with the 21 meV excitation lowering in energy 
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FIG. 4. (color online) Different AFM order used in the cal- 
culations. La atoms are light blue, O atoms green, As atoms 
purple and Fe atoms brown. The red and blue arrows show 
up and down spin, respectively, a) expanded non-magnetic 
unit cell, b) experimentally observed SDW, c) Checkerboard 
ordering d) striped ordering aligned ferromagnetically along 
the c-axis 



by approximately 11.9%. The intensity of the 32meV 
feature is half as strong as the nonmagnetic calculation. 

Fig. [2] compiles all of the experimental data and cal- 
culations of the different magnetic structures by showing 
several contour plots of the dynamical structure factor 
along (0, 0, L). Values of calculated intensities range from 



blue (low intensity) to red (high intensity) , and have been 
multiplied by the energy to improve visibility of the op- 
tical modes. The white dots show the experimentally 
determined frequencies with the intensity shown by the 
size of the dot. In each case, calculations with a mag- 
netic moment on the Fe show splitting between the two 
Ai g modes near 24 meV. At the zone center, the upper 
mode softens by 3.7 (SDW) - 5.3 (checkerboard) % and 
the lower mode softens by 9.1 (SDW) - 16.6 (checker- 
board) %. Calculated frequencies of these two modes 
are a few meV lower than observed. Comparing non- 
magnetic and spin-polarized calculations, the frequency 
of the zone boundary the upper La-As mode is essen- 
tially unchanged (< 0.6%). The intensity of the lower 
mode is strongest near (0,0,9), and the intensity of the 
upper mode is strongest near (0,0,8). The frequency 
of the lower mode differs by 2 meV. The SDW calcu- 
lation best matches the experimental frequency, but the 
checkerboard pattern best matches the observed splitting 
between these two modes. We note that the checker- 
board ordering also introduces a pronounced softening 
of the longitudinal acoustic mode when compared to the 
non-magnetic calculation and other magnetically ordered 
structures. Finally, we point out that the intensity of the 
optical mode near 10 meV is highest in nonmagnetic cal- 
culation. Overall changes in the phonon frequencies and 
intensities indicate the complex and subtle effects that 
magnetic ordering has on the lattice dynamics. 

Despite the changes introduced by magnetic order, all 
the spectra are qualitatively similar for different mag- 
netic orders, and in better agreement with experiment 
compared with nonmagnetic calculations. This might 
be understood to occur as a consequence of Fe mo- 
ments still being present above T/v, though without long- 
range orderP^ Compared to nonmagnetic calculations, 
imposing an AFM ordering better describes phonons in 
LaFeAsO. Consequently, it is likely that the presence of 
Fe moments, ordered or not, affects the force constants. 
Considering only z-polarized phonon branches contain- 
ing La and As motion significantly reduces the number of 
force constants that contribute. First, only the "zz" term 
in the 3x3 force constant tensor can contribute, greatly 
simplifying comparisons between different magnetic unit 
cells. Fe and O are essentially stationary in the modes 
considered, meaning force constants between Fe-Fe,Fe-0 
and 0-0 atoms do not contribute. La-La and La-0 force 
constants are essentially unchanged in each calculation, 
and the bond distance between La-Fe is large and the 
resulting force constant small. Therefore, we can limit 
ourselves to the "zz" term for La-As, As- As, and Fe-As 
force constants. Of these, the Fe-As force constant is 
the largest by an order of magnitude. Even with these 
simplifications, there was no clear softening of any spe- 
cific pair- wise force. This provides additional evidence 
for T. Yildirim's observation that changes in the phonon 
modes are due to a complicated renormalization rather 
than softening of a single pair-wise forceP In the non- 
magnetic calculation the on-site force constants, which 
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are a sum of all pair-wise interactions, are (in eV/A ) 
(11.2,11.2,8.9) for Fe and (10.5,10.5,9.6) for As. They 
show significant softening around 10-20% with the intro- 
duction of magnetic order, in good agreement with T. 
Yildirium's work. Small differences on the order of 0.2 
eV/A or less in the on-site force constants are likely 
from slightly different lattice parameters chosen in our 
calculations. 

In summary, we have measured the phonon dispersion 
along (0,0, L) tetragonal phase of LaFeAsO at room tem- 
perature, well above the magnetic ordering temperature 
of 138K. Nonmagnetic calculations fail to reproduce the 
observed splitting between two A\ g phonon modes at 22 
and 26 meV. Spin-polarized first-principles calculations 
imposing a number of hypothetical antiferromagnetic or- 
ders are qualitatively similar and in better agreement 
with the experimental results than non-spin-polarized 



calculations. The presence of Fe-spins are necessary to 
predict the observed spectrum above T/v, however the 
renormalization of the force constants is quite complex 
and cannot be reduced to a single pair-wise force con- 
stant. 
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